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I  INTRODUCTION 


Theoretical  work  in  wave  propagation  in  randomly  irregular  media 
has  progressed  considerably  over  the  last  several  years  [Flatte  et  al. , 
1980;  Rino  et  al. ,  1980b;  Rino,  1980]-  A  collection  of  comparatively 
simple  formulas  has  emerged  that  characterize  the  most  important 
channel  effects — temporal  coherence  loss,  frequency  coherence  loss,  and 
angle  scattering, — under  conditions  of  strong  scattering,  such  as  would 
occur  in  the  ionosphere  following  a  high-altitude  nuclear  detonation. 

To  the  extent  that  these  results  have  been  tested,  they  have  proven 
accurate . 

Insofar  as  basic,  theoretical  work  on  propagation  is  concerned 
little  more  is  required  for  predicting  system  performance  in  highly 
disturbed,  but  known  propagation  environments.  The  basic  problem  of 
accurately  characterizing  such  propagation  environments  remains, 
however,  and  new  issues  continue  to  emerge.  For  example,  in  the 
power-law  environments  that  seem  characteristic  of  all  turbulent 
channels,  the  severity  of  the  propagation  effects  are  sensitive  to 
changes  in  the  spectral  index.  To  study  such  subtle  effects,  a 
refined  propagation  theory  is  important  because  "channel  sounding" 
remains  the  most  economical  means  of  obtaining  data  on  irregularity 
structures.  Under  this  contract,  we  have  investigated  a  new  technique 
for  beacon  diagnostics. 

For  some  time  there  has  been  concern  about  the  accuracy  with 
which  the  scintillation  structure  "mirrors"  the  underlying  irregularity 
structure  that  causes  the  scintillation.  This  concern  was  intensified 
when  the  Wideband  satellite  data  showed  a  spectral  index  that  was 
not  only  shallower  than  was  expected,  but  also  variable.  This  dilemma 
has  only  recently  been  resolved  by  carefully  analyzing  data  from  the 
Atmospheric  Explorer  Satellite  E  (AE-E) [Livingston  et  al.,  1980]. 

These  data  showed  that  the  in-situ  spectral  index  is  indeed  more 


shallowly  sloped  than  expected  and  that  it  varies  systematically,  as 
the  Wideband  satellite  data  implied. 

In  addition  to  the  shape  of  the  spectral  density  function  (SDF), 
we  would  like  to  measure  the  absolute  electron-density  fluctuation 
level  or  turbulent  strength  (as  distinct  from  N/N)  for  each  spectral 
component.  To  do  so  from  conventional  phase-scintillation  data  involves 
a  two-fold  complication: 

(1)  The  level  of  phase-scintillation  depends  on  the  propagation 
geometry  relative  to  the  principal  irregularity  axis  of 
elongation  and  the  corresponding  "stretch"  factors  or 
axial  ration. 

(2)  The  level  of  phase-scintillation  depends  on  the  integral 
along  the  propagation  path. 

Because  the  irregularity  anisotropiiy ,  the  patli  length  of,  and,  in 
particular,  the  structure  variation  along  the  propagation  patli  are  a 
priori  unknown,  a  completely  unambiguous  determination  of  turbulent 
strength  from  such  scintillation  data  is  not  feasible. 

In  the  DNA  PLUMEX  experiment  conducted  at  Kwajalein  during  Julv- 
August  1979,  however,  a  downward-looking  radio  beacon  was  flown 
together  with  a  sophisticated  complement  of  in-situ  probes.  Unlike 
the  satellite-beacon  measurements,  in  the  rocket  configuration  the 
velocity  component  along  the  line  of  sight  is  much  larger  than  the 
transverse  component.  A  comparative  analysis  of  the  in-situ  probe 
and  beacon  data  from  PLUMEX  has  shown  that  there  is  a  simple  derivative 
relationship  between  the  beacon  phase  and  the  in-situ  density  that 
does  not  involve  any  unknown  parameters  [Rino  et  al . ,  1980;  Pet  r i ceks , 
1980].  Thus,  the  beacon  data  can  be  used  to  determine  the  turbulent 
strength . 

Because  of  the  comparative  simplicity  of  such  rocket-beacon 
measurements,  we  have  developed  a  tiieory  (Section  II  of  this  report) 


U 


1  1 


for  interpreting  rocket  beacon-phase  spectra  where  Vy  '•>  .  The 

results  show  that  as  long  as  3  =  vj^/v||  is  less  than  •  40  percent, 

4>z(k)  =  re2A2  $1(k)/k2  , 

where  4>  (k)  is  the  beacon-phase  SDF  and  <5>^(r)  is  the  one-dimensional 
SDF  as  measured  by  the  in-situ  probe.  The  relationship  obviously 
breaks  down  at  small  wave  numbers,  but  the  analysis  shows  that, 
depending  on  3,  the  relationship  holds  to  within  a  few  spectral 
resolution  cells  of  the  smallest  resolvable  Fourier  component. 

Because  the  classical  electron  radius,  r  ,  and  the  wavelength,  A, 
are  known,  the  beacon-phase  SDF  can  be  simply  and  unambiguously 
related  to  the  in-situ  SDF.  The  technique  is,  therefore,  potentially 
useful  for  future  diagnostic  rocket  probes. 

In  Section  III  of  this  report,  we  have  investigated  the  simplest 
form  of  the  mutual  coherence  function,  which  is  useful  for  modeling 
the  Doppler  spectrum  of  a  scintillating  signal.  An  approximate 
formula  is  developed  that  relates  the  perturbation  strength  to  the 
spectral  width  without  having  to  evaluate  gamma  functions.  Two 
limiting  analytic  forms  for  the  Doppler  spectrum  are  derived:  one 
characterizes  shallowly-sloped  phase  spectra  and  has  a  power-law  form; 
the  second  has  a  Gaussian  form  and  is  appropriate  for  more  steeply- 
sloped  spectral-density  functions.  Numerical  computations  are  presented 
that  show  the  variation  of  the  spectral  shape  with  changing  the 
spectral  index  between  the  two  extremes  where  analytic  results  can 
be  obtained. 

These  results  are  useful  for  both  predictive  modeling  and  data 
analysis.  The  Doppler  spectral  model  is  being  verified  concurrently 
using  Wideband  satellite  data. 
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II  RADIO-BEACON  MEASUREMENTS  OF  LOCAL  IRREGULARITY  STRUCTURES 

A.  Background 

Over  the  last  decade,  extensive  phase  scintillation  data  have 
been  obtained  from  low  orbiting  and  geosynchronous  satellites  [Crane, 
1977;  Davis  et  al.,  1975;  Rino  et  al . ,  1980b].  In  such  measurements 
both  the  signal  source  and  receiver  are  generally  well  removed  from 
the  disturbed  medium.  Ignoring  diffraction  effects  and  temporal 
changes  in  the  medium,  the  instantaneous  phase  perturbation  is  given 
by  the  formula 


£ 

P 

<p(t)  =  kR(t)  -  rgA  J’  N  /v^t ,  n 

0  ' 

where  k  =  2tt/A,  R(t)  is  the  range  to  the  satellite,  1  is  the  length 
of  the  propagation  path  within  the  disturbed  medium,  and  Ne(p,  z)  is 
the  local  electron  density. 

The  path  length  changes  with  time,  but  l  is  assumed  nearly 

P 

constant  over  periods  long  compared  to  the  time  scale  for  phase 
changes  of  interest.  If  v  =  v  a  ,  the  phase  SDF  takes  the  form 

i  y  Y 


idn 


(1) 


<Kf) 


G 


where  '{’(k,  <z)  is  the  three-dimensional  SDF  of  Ng(p,  z).  The 
corresponding  one-dimensional  SDF  measured  by  an  in-situ  probe  is 


f  C  / 


\  dK'  dK 


h* 
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It  is  readily  shown  [Rino,  1979]  that  if  <!>(?,  k^)  q  (^v+l)^  then 
( f )  ^  f  and  4>^(f)  “  f  ^  .  Thus,  if  v  =  1.5,  then  the  phase 

_3 

SDF  has  the  power-law  form  f  ,  whereas  the  one-dimensional  in-situ 

-2 

has  the  power-law  form  f  .  The  difference  of  unity  between  the  in-situ 
and  phase  spectral  indices  has  been  verified  by  comparing  spectral 
measurements  derived  from  Wideband  satellite  phase  data  and  AE-E 
in-situ  data  [Livingston  et  al.,  1980]. 

For  a  downward-looking,  rocket-borne  radio  beacon,  the  signal 
phase  takes  the  form 


<J>(t)  = 


k(z0  +  ''z0 


r  A 
e 


/ 


z  +v  t 
0  z 


Ne(v^  t,n)dn 


(A) 


where  v^  and  v^  are  the  relative  velocity  components  along  and 
transverse  to  the  line  of  sight  to  the  receiver.  The  fixed  path 
length  Zq  effectively  establishes  the  initial  condition  for  the 
indefinite  integral  that  vzt  maps  out. 

In  the  early  development  of  guided  missiles,  radio-beacon 
measurements  were  used  for  trajectory  evaluation.  As  with  satellite 
radio  navigation,  the  ionospheric  contribution  was  a  source  of  error. 
As  the  techniques  were  refined,  however,  the  radio-beacon  data  came  to 
be  used  to  determine  ionospheric  density  profiles  [ Berning,  1951]. 

If  v  is  small,  we  have  only  to  remove  the  geometric  Doppler  term 
kv  t  in  Eq.  (A)  and  differentiate  the  residual  to  obtain  N  (z). 

Note  that  the  measurement  gives  absolute  electron  density  because 
only  phase  changes  are  involved. 


More  recently,  the  technique  has  been  exploited  to  measure 
small-scale  irregularity  structure.  In  July  1979,  a  rocket-borne 
radio  beacon  was  launched  from  Roi  Namur  Island  in  the  Kwajalein 
Atoll  into  a  highly  disturbed  equatorial  ionosphere.  The  results 
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of  the  beacon  and  probe  structure  comparisons  are  described  in 
Rino  et  al .  [1980c].  The  data  interpretation  was  based  on  the 
assumption  that  in  Eq.  (4)  could  be  disregarded.  This  assumption 
was  justified  ex  post  facto  by  the  agreement  between  the  in-situ  and 
beacon  data. 

In  this  section  we  develop  the  theory,  in  detail,  to  assess 
the  general  applicability  of  the  radio-beacon  measurement  technique 
for  determining  local  in-situ  irregularity  structures.  We  begin 
by  computing  the  general  relation  between  the  beacon-phase  SDF  and 


B.  Spectral  Relations  for  Rocket-Beacon  Phase  Measurements 

To  simplify  notation,  we  let  7.  =  v^t  and  consider  the  indefinite 
integral 


Nt(z) 


■/ 


Vz 


Ne((-.z, 


n)dn 


(5) 


where 


(■'  =  v 

1 


/v 

z 


(6) 


If  we  assume  that  NO,  7.)  can  be  modeled  by  a  statistically  homogeneous 
random  process  throughout  the  measurement  region,  we  can  write 


N  (.  ,  7.)  = 


■V' 


(7) 


where 


<dV  •  •'»  -  —  (•*.  -,)«-•  -  -  >z)  <»> 


is  a  purely  formal  definition  of  the  orthogonal  increments  property  of 

the  Fourier  spectrum  d£_,(K,  >  ). 

N  z 

By  substituting  Eq.  (7)  into  Eq .  (5)  and  changing  the  order  of 
integration,  we  can  eliminate  the  indefinite  integral  in  Eq.  (5). 

The  result  is 


Vz)  =  jj j  oxp{-h 


exp{- i k  (z  +  z0> >-l 


Because  we  can  only  compute  the  finite  Fourier  transform,  we  shall  only 
consider  the  integral 


V  ?)  =f  NT(z)exp{i.  7z^dz 


In  fact,  N^(z)  must  be  modified  to  avoid  contamination  of  the  spectral 
estimate  because  of  "end  point  mismatch."  Such  subtleties  need  not  be 
formally  carried  through  the  analysis,  however. 

If  Eq.  (9)  is  now  substituted  into  Eq .  (10)  and  the  integration 
over  z  performed,  the  result  is 


W  f f /™pC-«W  v  z0)  — N— 


9 


where 


W  q2’  q3;  V  = 


exp<-i(q1  +  q2  -  q3>L /  -  1  /  j 

(q1+q2-q3) - 


exp 


{-Kqj^  -  q3)L}  ~  1 


(q3  -  q3> 


(12) 


An  estimate  of  the  SDF  of  NT(z)  is  obtained  by  computing 
|NT«z)|2/L  and  averaging  to  smooth  the  statistical  fluctuations. 
The  result  is  close  to  the  expectation  value. 


!z(k7)  =  <|nt(kz)|Z>/L 


(13) 


To  evaluate  Eq .  (13),  we  substitute  NT(<z)  from  Eq .  (11)  and 
use  the  orthogonal  increments  property  Eq .  (8)  to  obtain 


(k 

z  z 


r/> 

>  v  Zq)I 


2  i(K;  k?  dK<  ^z 

k'2  (2-i) 2  2n 

z 


(14) 


We  note  that  there  is  no  singularity  at  k'  =  0  because 


1  im 

V°  V*'**’  V  Kz;  Z0)kz 


_2  s in2 [ (if  '  •  3  -  kz)L/2] 
(k '-3  -  k  )2 


(15) 


Let  us  first  consider  the  special  case  for  which  8  =  0.  Then, 
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(16) 


,  12  VKz}  dKz 

:(lCz)  j  L  |DL(0,  V  V  VI  y.2  2ti 


where 


CO 


dK 

(2 7T  )2 


(17) 


is  the  one-dimensional  spatial  SDF  derived  from  a  probe  scanning 
along  the  z  axis  which  should  not  be  confused  with  the  temporal  SDF 
defined  by  Eq.  (3). 


The  general  behavior  of  Eq.  (16)  is  readily  established.  We 

first  isolate  the  k'  regime  where  Eq.  (15)  applies.  It  is  sufficient 
z 

to  take  k'  <  2tt/L,  whereby 


(18) 


This  contribution  to  Eq.  (16)  becomes  vanishingly  small  as  ^  increases; 
however,  for  large  k 
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1 

L 


V°’ 


?  *  2 
^  sin 

(k  -  k') 

- 

.  .  r,  \ 

.  z  z  2] 

<*z  - 


09) 


Now,  as  long  as  the  variation  of  i>  (r ')  is  small  over  k '  intervals, 

1  z  z 

comparable  to  2n/L,  we  can  show  that 


(20) 


[ Rosenblatt ,  pp.  169-180,  1962],  which  is  sometimes  referred  to  as  the 

2 

"sifting  property"  of  sine  . 

As  long  as  k  ■  2 r: / L ,  the  =  0  behavior  of  Eq .  (14)  is 
summarized  by 


1  (r  )  =  0.  )hA 

z  z  1  z  z 


(21) 


The  relationship  Eq.  (21)  is  useful  because  it  involves  no  unknown 
paramters,  e.g.,  the  length  of  the  propagation  path  or  the  axial 
ratios,  which  characterize  the  anisotropy  of  the  medium. 

At  the  opposite  extreme  where  fl  is  very  large,  we  can  disregard 
q2  in  Eq.  (12)  where  is  appears  next  to  q^.  The  result  is 
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i.i ip'ipwi in 4  ^pi i 


-v 

In  the  usual  approximation  the  integral  over  k  is  replaced  by  z  1 (t  ;  0) . 

*  2  U  2 
Finally,  by  letting  (-;  =  pa^  and  applying  the  sifting  property  of  sine 

to  the  k  integration,  we  have 


0) 


(23) 


which  is  essentially  equivalent  to  Eq .  (2).  Note  that  z in  Eq .  (23) 
is  the  distance  from  the  point  where  the  disturbance  starts  and  is, 
therefore,  unknown. 

For  0  P  1,  which  is  the  case  of  primary  interest,  we  first 
isolate  the  region  near  =  0  and  use  Eq.  (15).  In  place  of  Eq.  (18) 
we  now  have 
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>  2 
If  we  let  B  =  Ba^  and  using  the  sifting  property  of  sine  ,  the  right- 

hand  side  of  Eq.  (24)  can  be  replaced  by 


(25) 


For  a  power-law  SDF  this  contribution  again  decreases  rapidly  with 
increasing  k 

The  behavior  of  Eq.  (14)  for  the  region  outside  |k'|<2h/L  can  be 
determined  with  the  aid  of  Figure  1  where  the  k'  -  k'  plane  has  been 

y  Z 

partitioned  into  two  regions  separated  by  the  line  k'  =  k'/B.  To  the 

y  z 

right  of  this  line. 


sin' 


6  +  k;  -  kz)L/2 I 


(k'-|  +  k'  -  K  )2L 
z  z 


(26) 


14 


T 


j  mi  4wu- 


In  Figure  3,  the  region  where  Eq.  (26)  breaks  down  is  evident.  The 
effect  of  increasing  3  is  shown  in  Figure  4. 

We  have  already  established  that  the  contribution  to  Eq.  (14) 
foi:  small  k'  is  negligible  when  <<  2r/L.  Thus,  there  is  a  minimum 

value  of  k  as  shown  in  Figure  1.  Beyond  an  arc  drawn  from  the 
z 

intersection  of  k'  =  2tr/L  and  =  (k  -  k')/3,  the  contribution  to 
z  y  min 

Eq.  (14)  of  <t> (k ^  k')/k"  must  be  negligible.  We  can  then  apply  the 
sifting  property  of  Eq.  (26)  to  obtain 


The  second  approximation  applied  because  k  <<  k  '  •  £  in  the  region 
where  Eq.  (26)  is  valid.  For  a  typical  power-law  SDF ,  the  validity 
conditions  for  Eq.  (27)  are  easily  achieved  3  <-  *  4. 
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Ill  DOPPLER  SPECTRUM 


A.  Background 

For  some  predictive  scintillation  codes,  it  is  useful  to  model 
the  Doppler  spectrum  of  the  complex  signal,  v(t).  The  Doppler  spectrum 
is  formally  the  Fourier  transform  of  the  mutual  coherence  function 


R  (t)  =  <v(t)v  (t  +  x)^ 


(28) 


Under  reasonable  assumptions,  moreover,  R  (t)  admits  the  simple 
representation 


Rv(t>  =  exp{-  \  D,S|  (veff  T)|  ,  (29) 

where 

Ds  ,  (y)  =  2n(V2  [l  -  (.A;  (y)]  (30) 

is  the  phase  structure  function,  o  is  the  rms  phase,  and  .  ,  (y) 
is  the  normalized  phase  autocorrelation  function  [i^.(0)  =  1]. 


Equation  (29)  is  strictly  valid  only  in  the  plane  normal  to  the 
pronagation  direction,  that  is,  when  parallel  propagation  effects 
are  negligible.  The  correction  factor  that  must  be  applied  to 
accommodate  parallel  propagation  effects  are  described  in  Rino, 
Section  V,  [1930];  thus,  we  shall  only  consider  Eq.  (29)  here. 

If  the  three-d imens ional  spatial  SDF  of  the  irregularities  lias 
the  general  power-law  form 


abC 
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can 


then  ,  s.(y),  the  two-dimensional  Fourier  transform  of  Eq.  (31), 
be  evaluated  as 


,  .  (y) 


2 

(  -  1/2) 


-1/2 


K  -l/2(q0v) 


(32) 


where  y  is  a  quadratic  form  in  the  spatial  separation  coordinates, 
and  .  ,  as  described  in  Rino  |  1979a,  b|.  However,  in  spite-  of 

X  V  “ 

the  comparatively  simple  torm  of  Eq .  (32),  it  is  too  cumbersome  for 
direct  use  in  Eq.  (29),  and  we  seek  simpler  approximate  forms. 

The  simplest  approximation  that  retains  a  dependence  on  the 
spectral  index  parameter,  ,  is 


>  .  9  _  ] 

n  (v)  sc./.v'“ 

2 

where  C  is  the  phase  structure  constant 


(33) 


C, 


2 


2  2,,  5^ _ 2.  (1. 3  -  v) _ 

°  P  2  (0.3  +  )(2v  -  1)22  _1 


(34) 


which  is  valid  for  j  1.5.  Equation  (33)  is  obtained  as  an 
asymptotic  approximation  Eq.  (30)  in  the  limit  as  q^  becomes  very 
small  [ R ino ,  1979b].  Note  that  Eq.  (33)  does  not  retain  an  explicit 
dependence  on  q^. 

For  more  steeply  sloped  spectra  (  1.5),  a  quadratic 

approximation  of  the  form 


I) 


(35) 
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is  admissible.  In  Eq.  (35),  is  the  first  nonzero  coefficient  in 
a  Taylor  series  expansion  of  D^(y).  The  expansion  coefficients 
generally  do  not  exist-  independently  from  an  inner-scale  cut-off, 
f Rino  et  al.,  1980],  although  is  independent  of  q^  and  for 
v  >  1.5. 


For  the  entire  range  of  v  values,  R^(x)  admits  the  approximate 


form 


where 


R,(T)  -  expj-  1/2  C^2  |yp»<2v-l,2)j 


(36) 


-zqj  -  y± 


r(0.5  +  v)(2v  -  l)v 


2v-l 


0.5  <  v  <  1.5 


2nC 


log(qi/q() 


v  =  1.5 


(37) 


r(v  -  1.5) 
,4T(v  -  0.5) 


v  >  1.5 


2 

A  plot  of  2rrn  ,  /C  as  a  function  of  v  is  illustrated  in  Figure  5. 

0(f)  p 

The  singular  behavior  at  v  =  1.5  is  indicative  of  a  breakdown  in  both 
the  asymptotic  and  quadratic  approximations.  The  dashed  line  is 
probably  closer  to  the  correct  functional  dependence  of  an  asymptotic 
approximation  of  the  form  of  Eq.  (32). 
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For  predict  ivt?  modeling,  therefore,  it  is  convenient  to  Laki 


2  C  2/C  S. 2  -4. 5.  +5. 5  1.0  <  <2.0 

P 


which  follows  the  dashed  curve  in  Figure  3  very  closely.  Approximations 

sucii  as  Eq.  (38),  arc  part icularly  convenient  because  of  the  evidence 

that  the  spectral  index  itself  varies  systematically  with  changing 

turbulent  strengtli  [Livingston  et  al.,  1080:,  Rino  et  a].,  1980c], 

In  any  case,  for  our  purposes  here,  we  are  only  interested  in  the 

behavior  of  Eq .  (36)  as  a  function  of  .  Because  of  the  simple  power- 

2 

law  form  of  the  argument  of  the  exponential  function,  C  is 
effectively  only  a  scale  factor. 


B.  Special  Cases  and  Numerical  Computations 

To  compute  the  Doppler  spectrum,  we  must  evaluate  the  integral 


F  (f )  Ry(')  exp  |-2  ifijdx 


Substituting  y  =  v  ^r  into  Eq.  (36)  gives 


.min (2 


Vr> 


where 


0  vt.ff  l  c  2  /  min (2  -1,2) 


It  -  1,  Eq .  (39)  admits  the  exact  solution 

:(f)  =1 - 2 -  .  (42) 

0  1  +  (2/n0f)Z 

Similar  l  v,  il  =  1.3,  tlien  1(f)  takes  the  Gaussian  form 

(f )  =  exp 

0 

Thus,  the  genera  1  trend  is  from  shallow  power-law  to  Gaussian  as 
ine reuses . 

This  is  verified  in  the  numerical  evaluations  of  Eq.  (39)  using 
the  fast  Fourier  transform  algorithm.  In  Figure  6,  t^T(f)  is  plotted 
for  =0.8,  =0.9,  and  =  1.0.  For  such  shallowly  sloped  phase 

spectra,  the  overall  shape  of  ! (f)  changes  little.  As  the  slope 
increases,  however,  shape  changes  rapidly — as  shown  in  Figure  7  where 
.^•‘  (f)  is  plotted  for  v  =1.3,  .'  =  1.4,  and  v  =  1.5. 

For  predictive  modeling,  therefore,  we  can  use  Eq .  (42),  where 

the  spectral  index  is  comparatively  small,  and  Eq .  (43)  for  the 

opposite  extreme.  It  should  he  kept  in  mind,  however,  that  1q  itself 
depends  on  the  spectral  index  parameter,  v ,  as  well  as  on  the 
perturbation  strength. 


,f) 


(43) 
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FIGURE  7  DOPPLER  SPECTRUM  FOR  MORE  STEEPLY  SLOPED 
SPECTRA 


2 


IV  CONCLUSIONS 


In  this  report  we  developed  a  theory  to  interpret  beacon-phase 
data  when  the  line-of-sight  velocity  component  is  much  larger  than  the 
transverse  component.  In  that  situation,  there  is  a  simple  relationship 
between  the  phase  spectral-density  function  and  the  one-dimensional 
in-situ  spectral-density  function  that  a  rocket-borne  density  probe 
would  be  measured.  Thus,  noth  the  spectral  shape  and  the  turbulent 
strength  can  be  determined. 

Because  the  rocket  beacon  signal  itself  is  the  data  source,  no 
telemetry  is  required.  Thus,  such  beacon  measurements  provide  a 
comparat ively  inexpensive  means  for  obtaining  quantative  ionospheric 
structure  measurements. 

Also,  a  simple  mathematical  model  was  developed  for  predicting 
the  Doppler  spectrum  of  a  satellite  signal  received  in  the  usual 
situation  in  which  the  transverse  Doppler  component  dominates.  This 
model  is  useful  both  for  data  interpretation  and  system  analysis. 

Host  other  theoretical  results  of  interest  have  been  verified. 

The  mutual  coherence  function  model  itself,  which  is  the  starting 
point  for  the  Doppler  spectrum  computations,  has  not  yet  been  tested 
against  real  data,  but  this  effort  is  being  pursued  under  a  separate 
contract . 
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ATTN:  MS  670,  J.  Hopkins 

Sandia  National  Laboratories 

Livermore  Laboratory 

ATTN:  B.  Murphey 
ATTN:  T.  Cook 

Sandia  National  Lab 

ATTN:  ORG  1250,  W.  Brown 
ATTN:  ORG  4241,  T.  Wright 
ATTN:  Space  Project  Div 
ATTN:  3141 
ATTN:  D.  Thornbrough 
ATTN:  D.  Dahlgren 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Aerospace  Corp 

ATTN:  I.  Garfunkel 
ATTN:  S.  Bower 
ATTN:  D.  Olsen 
ATTN:  N.  Stockwell 
ATTN:  R.  Slaughter 
ATTN:  V.  Josephson 
ATTN:  T.  Salmi 
ATTN:  1.  Straus 

University  of  Alaska 
ATTN:  T.  Davis 
ATTN:  N.  Brown 
ATTN:  Technical  Library 

Analytical  Systems  Engineering  Corp 
ATTN:  Radio  Sciences 

Analytical  Systems  Engineering  Corp 
ATTN:  Security 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 
Barry  Research  Corporation 


General  Electric  Co 

ATTN:  A.  Steinmayer 
ATTN:  C.  Zierdt 

General  Electric  Co 

ATTN:  F.  Reibert 

General  Electric  Tech  Services  Co,  Inc 
ATTN:  G.  Mi  liman 


ATTN:  E.  Knick 

Horizons  Technology,  Inc 
ATTN:  R.  Kruger 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 
HSS,  Inc 


ATTN:  J.  McLaughlin 

ATTN:  D.  Hansen 

BDM  Corp 

IBM  Corp 

ATTN:  L.  Jacobs 

ATTN:  F.  Ricci 

ATTN:  T.  Neighbors 

Berkeley  Research  Associates,  Inc 

University  of  Illinois 

ATTN:  Security  Supervisor  fo 

ATTN:  J.  Workman 

Betac 

Institute  for  Defense  Analyses 

ATTN:  J.  Bengston 

ATTN:  J.  Hirsch 

ATTN:  E.  Bauer 

Boeing  Co 

ATTN:  H.  Wo If hard 

ATTN:  J.  Aein 

ATTN:  M/S  42-33,  J.  Kennedy 

ATTN:  G.  Hall 

International  Tel  &  Telegraph  Corp 

ATTN:  S.  Tashird 

ATTN:  G.  Wetmore 

ATTN:  Technical  Library 

Booz-Allen  &  Hamilton,  Inc 

ATTN:  B.  Wilkinson 

JAYCOR 

ATTN:  J.  Sperling 

University  of  California  at  San  Diego 

ATTN:  H.  Booker 

JAYCOR 

ATTN:  J.  Doncarlos 

Charles  Stark  Draper  Lab,  Inc 

ATTN:  D.  Cox 

Johns  Hopkins  University 

ATTN:  J.  Gilmore 

ATTN:  T.  Potemra 

Communications  Satellite  Corp 

ATTN:  J.  Phillips 

ATTN:  T.  Evans 

ATTN:  D.  Fang 

ATTN:  J.  Newland 

ATTN:  P.  Komiske 

COMSAT  Labs 

ATTN:  R.  Taur 

Kaman  Tempo 

ATTN:  G.  Hyde 

ATTN:  DAS I AC 

Cornell  University 

ATTN:  W.  McNamara 

ATTN:  T.  Stephens 

ATTN:  M.  Kelly 

ATTN:  W.  Knapp 

ATTN:  D.  Farley,  Jr 

Electrospace  Systems,  Inc 

Linkabit  Corp 

ATTN:  I.  Jacobs 

ATTN:  H.  Logston 

ESL,  Inc 

Litton  Systems,  Inc 

ATTN:  R.  Grasty 

ATTN:  J,  Marshall 

General  Electric  Co 

Lockheed  Missiles  A  Space  Co,  Inc 
ATTN:  W.  Imhof 

ATTN:  M.  Bortner 

ATTN:  M.  Walt 

ATTN:  A.  Harcar 

ATTN:  R.  Johnson 

Lockheed  Missiles  A  Space  Co,  Inc 
ATTN:  Dept  60-12 

M.  I  .T.  Lincoln  Lab 

ATTN:  D.  Towle 

Martin  Marietta  Corp 

ATTN:  R.  Heffner 

McDonnell  Douglas  Corp 


General  Research  Corp 

ATTN: 

G. 

Mroz 

ATTN:  J.  Ise,  Jr 

ATTN: 

N. 

Harris 

ATTN:  J.  Garbarino 

ATTN: 

W. 

Olson 

ATTN: 

J. 

Moule 

Harris  Corp 

ATTN: 

R. 

Halprin 

Meteor  Conmuni cat ions  Consultants 
ATTN:  R.  Leader 


DEPARTMENT 

OF  DEFENSE  CONTRACTORS 

(Continued! 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Mission  Reserach  Corp 

Rockwell  International  Corp 

ATTN 

R.  Ki lb 

ATTN 

R.  Buckner 

ATTN 

R.  Bogusch 

ATTN 

0.  Sappenfield 

Rockwell  International  Corp 

ATTN 

F.  Fajen 

ATTN 

S.  Quilici 

ATTN 

R.  Hendrick 

ATTN 

S.  Gutsche 

Santa  Fe  Corp 

ATTN 

Tech  Library 

ATTN 

D.  Paolucci 

Mitre  Corp 

Science  Applications,  Inc 

ATTN 

G.  Harding 

ATTN 

C.  Smith 

ATTN 

A.  Kynriel 

ATTN 

D.  Hamlin 

ATTN 

B.  Adams 

ATTN 

E.  Straker 

ATTN 

C.  Callahan 

ATTN 

L.  Linson 

Mitre  Corp 

Science  Applications,  Inc 

ATTN 

M.  Horrocks 

ATTN 

SZ 

ATTN 

J.  Wheeler 

ATTN 

W.  Foster 

Science  Applications,  Inc 

ATTN 

W.  Hall 

ATTN 

J .  Cockayne 

Pacific-Sierra  Research  Corp 

SRI  International 

ATTN 

E.  Field,  Jr 

ATTN 

W.  Jaye 

ATTN 

F.  Thomas 

ATTN 

R.  Leadabrand 

ATTN 

H.  Brode 

ATTN 

D.  Neil  son 

ATTN 

J.  Petrickes 

Pennsylvania  State  University 

ATTN 

W.  Chesnut 

ATTN 

Ionospheric  Research 

Lab 

ATTN 

R.  Livingston 

ATTN 

R.  Tsunoda 

Photometries,  Inc 

ATTN 

G.  Smith 

ATTN:  I.  Kofsky 

ATTN 

G.  Price 

ATTN 

M.  Baron 

Physical  Dynamics,  Inc 

ATTN 

A.  Burns 

ATTN 

E.  Fremouw 

4  cy  ATTN 

C.  Rino 

Physical  Research,  Inc 

Syl vania  Systems  Group 

ATTN 

R.  Deliberis 

ATTN 

M.  Cross 

RAD  Associates 

Technology  International  Corp 

ATTN 

R.  Turco 

ATTN 

W.  Boquist 

ATTN 

C.  Greifinger 

ATTN 

M.  Gantsweg 

TRI-COM,  Inc 

ATTN 

B.  Gabbard 

ATTN 

D.  Murray 

ATTN 

W.  Wright 

ATTN 

R.  Lelevier 

TRW  Defense  &  Space  Sys  Group 

ATTN 

F.  Gilmore 

ATTN 

D.  Dee 

ATTN 

W.  Karzas 

ATTN 

R.  Plebuch 

ATTN 

H,  Ory 

ATTN 

P.  Haas 

Utah  State 

University 

ATTN 

L.  Jensen 

R  &  D  Associates 

ATTN 

J.  Dupnik 

ATTN 

B.  Yoon 

ATTN 

K.  Baker 

Rand  Corp 

Visidyne,  Inc 

ATTN 

C.  Crain 

ATTN 

C.  Humphrey 

ATTN 

E.  Bedrozian 

ATTN 

J.  Carpenter 

Riverside  Research  Institute 

ATTN 

V.  Trapani 

35 


